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We show that phase separation can occur in a one-component liquid outside its coexistence curve 
(CX) with addition of a small amount of a solute. The solute concentration at the transition 
decreases with increasing the difference of the solvation chemical potential between liquid and gas. 
As a typical bubble-forming solute, we consider O 2 in ambient liquid water, which exhibits mild 
hydrophobicity and its critical temperature is lower than that of water. Such a solute can be expelled 
from the liquid to form gaseous domains while the surrounding liquid pressure is higher than the 
saturated vapor pressure pcx. This solute-induced bubble formation is a hrst-order transition in 
bulk and on a partially dried wall, while a gas film grows continuously on a completely dried wall. 
We set up a bubble free energy AG for bulk and surface bubbles with a small volume fraction (f>. It 
becomes a function of the bubble radius R under the Laplace pressure balance. Then, for sufficiently 
large solute densities above a threshold, AG exhibits a local maximum at a critical radius and a 
minimum at an equilibrium radius. We also examine solute-induced nucleation taking place outside 
CX, where bubbles larger than the critical radius grow until attainment of equilibrium. 

PACS numbers: 64.75.Cd Phase equilibria of fluid mixtures, including gases, hydrates, etc. 

82.60.Nh Thermodynamics of nucleation 

51.30.-|-i Thermodynamic properties, equations of state 


I. INTRODUCTION 

Recently, much attention has been paid to the forma¬ 
tion of small bubbles, sometimes called nanobubbles, in 
water They have been observed with a dissolved 

gas on hydrophobic surfaces and in bulk [l3l - [l8j| in 

ambient conditions (around 300 K and 1 atm), where the 
pressure in the bulk liquid region is larger than the satu¬ 
rated vapor pressure Pcx or outside the coexisting curve 
(CX). Their typical radius R is of order 10 — 100 nm and 
their life time is very long. The interior pressure is given 
by 2a(R ~30 atm for a bubble with R = 50 nm from 
the Laplace law, where a is the surface tension equal to 
72 erg/cm^. Strong attractive forces have also been mea¬ 
sured between hydrophobic walls in water due to bubble 
bridging (^-0,0, HI] . These effects are important in var¬ 
ious applications, but the underlying physics has not yet 
been well understood. 

In this paper, we ascribe the origin of bubble formation 
to a hydrophobic interaction between water and solute 
[l9l - l^ . In our theory, the solute-induced phase separa¬ 
tion generally occurs in equilibrium when the solvent is 
in a liquid state outside CX and the solute-solvent in¬ 
teraction is repulsive. Most crucial in our theory is the 
solvation chemical potential of the solute /is(n,T) de¬ 
pending on the solvent density n and the temperature 
T. With increasing such a repulsion, its difference A/Hs 
between the liquid and gas phases can be considerably 
larger than the thermal energy ksT (per solute parti¬ 
cle). In this situation, the solute molecules are repelled 
from the liquid to form domains of a new phase (in gas, 
liquid, or solid). Supposing bubbles with a small vol¬ 
ume fraction </>, we set up a free energy AG accounting 
for considerably large A/Xs/fcsT'. Then, its minimization 
with respect to (j) and the interior solute density n[ yields 


the equilibrium conditions of bubbles in liquid (those of 
chemical equilibrium and pressure balance). 

As a bubble-forming solute in water, we treat O2, 
which is mildly hydrophobic with AfXs/ksT = 3.44 on 
CX at T = 300 K. Furthermore, the critical tempera¬ 
ture and pressure of water and O2 are given by (647.3 
K, 22.12 MPa) and (154.6 K, 5.043 MPa), respectively. 
Notice that the critical temperature of water is much 
higher than that of O2 (and than those of N2, H2, and 
Ar etc) due to the hydrogen bonding in water. As a 
result, no gas-liquid phase transition takes place within 
bubbles composed mostly of O2 in liquid water in ambi¬ 
ent conditions. In contrast, strongly hydrophobic solutes 
usually form solid aggregates in liquid water except for 
very small solute densities [Tol - I^ . 

In our theory, solute-induced bubbles can appear out¬ 
side CX only when the solute density exceeds a thresh¬ 
old density, where the threshold tends to zero as the liq¬ 
uid pressure approaches Pcx- In particular, above the 
threshold density, a surface bubble (a gas film) appears 
on a hydrophobic wall in the temperature range T < Td 
(T > Td). As is well known, this is possible only on CX 
without solute (in one-component fluids). Here, Td is the 
drying temperature [2^, [2^ determined by the solvent- 
wall interaction, so it is insensitive to a small amount of 
solute with mild A/ig. With a solute below the threshold 
density, we predict only a microsopically depleted layer 
outside CX (as in one-component fluids). Indeed, some 
groups [23-I23 detected only microscopic depletion layers 
on a hydrophobic wall with a dissolved gas, while other 
groups observed surface bubbles [Ml. 

On the other hand, to prepare stable bulk bubbles, 
macroscopic gas bubbles composed of O2 etc, have been 
fragmented by stirring in liquid water |15l4l8l| . In such 
measurements, Ohgaki et al.[ll realized bubbles with 
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i? ~ 50 nm in quasi-steady states, where the bubble den¬ 
sity was Ub ^ 19 /im“^ and the bubble volume fraction 
was (j) ~ 0.01. We shall see that the nucleation barrier of 
creating solute-induced bulk bubbles in quiescent states 
is too high for nucleation experiments for a gas such as 
O 2 except for very high liquid pressures. 

As a similar bulk phenomenon, long-lived hetero¬ 
geneities have also been observed in one-phase states of 
aqueous mixtures with addition of a salt or a hydropho¬ 
bic solute [sol - ls^ . Dynamic light scattering experiments 
indicated that their typical size is of order 100 nm. Theo¬ 
retically, such a phase separation can occur if the solute- 
solvent interaction is highly preferential between the two 
solvent components [3ll. 

This paper is organized as follows. In Sec.II, we will 
present a thermodynamic theory of bubbles induced by 
a small amount of solute, where the liquid pressure and 
the total solvent and solute numbers are fixed. In Sec.Ill, 
we will set up a bubble free energy AG. In Sec.IV, we 
will examine solute-induced nucleation. In addition, in 
Sec.me and Appendix A, we will briefly examine bubble 
formation at fixed chemical potentials and at fixed cell 
volume. 


II. EQUILIBRIUM BUBBLES WITH 
HYDROPHOBIC SOLUTE 

We consider a one-component solvent, called water, in 
a liquid state outside the coexistence curve (CX). We 
then add a small amount of a neutral, hydrophobic so¬ 
lute (impurities). The total solvent and solute numbers 
are fixed at N = Vn and Ni = Vfii, respectively, with fi 
and hi being the initial water and solute densities. Here, 
h is larger than the liquid density on CX before bub¬ 
ble formation. We keep the pressure in the liquid region 
at the initial value p larger than the saturated vapor (co¬ 
existence) pressure Pcx by attaching a pressure valve to 
the cell, as illustrated in Fig.l. We do not assume the 
presence of surfactants and ions (see remarks in Sec.IVA). 


A. Solvation chemical potential and Henry’s law 

We assume that the molecular volume of solute vi is 
of the same order as that of solvent Uw, since large hy- 
drophobic impurities tend to form solid precipitates |22l - 
[ 2 ^. We then consider the Helmholtz free energy density 
/ depending on the water density n and the solute den¬ 
sity ni in the dilute limit of solute. In this paper, we 
neglect the solute-solute interaction to obtain 

/(n,ni) = /w(n) -|- kBTni[\n{niVi) - 1 -I- i^s(n)], (1) 

where fwin) is the Helmholtz free energy density of pure 
water and iZs(n) is related to the solvation chemical po¬ 
tential Psin) in the limit of small ni by 


( 2 ) 


Hereafter, the T-dependence of the physical quantities 
will not be written explicitly. 

Note that the combination ln[?;i/Aj] -I- z/g can be de¬ 
termined unambiguously in thermodynamics in the limit 
ni —S' 0, where Ai is the thermal de Broglie length 
(cx T' -1/2). Thus, vi may be chosen to be independent 
of n without loss of generality. It is known that the en- 
tropic contribution to Vg is crucial for nonpolar impurities 
in liquid water (2ll - l^ . 

From eq.(l) we calculate the chemical potential of wa¬ 
ter p and that of solute pi as 


k- = dfyj/dn + ksTnigsin), 

Pi = kBT[\n{nivi) + i^s(n)], 

where we define 

dvg 1 5^1 

~ dn~ kBT dn ■ 

The pressure p = np + nipi — / is written as 


( 3 ) 

( 4 ) 

( 5 ) 


p = [ndf^/dn - f^] -|- fcsTni)! -I- ngs{n)], (6) 


where the first term is the contribution from the solvent 
and the second term from the solute. The typical size of 
Psin) is of the order of the solute molecular volume wj. 
In the presence of bubbles, p and pi take common values 
in gas and liquid, while the pressure in the bubbles is 
higher than that in the liquid by 2a jR. 

First, the homogeneity of pi in equilibrium yields 

ni = Uj exp[-i/g(n)], (7) 

as a function of n in two-phase states, where nj = 
exp{pi/kBT)/vi is a constant. In gas-liquid coexistence, 
let the water and solute densities be n' and n[ in gas and 
be h and fii in liquid, respectively. Then, eq. (7) gives 

hi/n( = exp[-Az/s], (8) 

where the Avg = Vs{n) — Vs{n'). This density ratio is 
called the Ostwald coefficient, which represents solubility 
of a gas [TqI - I^ . It is much smaller than unity for large 
Avg. Near CX, we may approximate Avg by its value on 
CX expressed as 

Avg = vg{ni^)-vg{nl^), (9) 

where the liquid and gas densities on CX 

of pure water. 

It is worth noting that Avg in eq.(9) is related to the 
Henry constant fcn [H, • From partitioning of a solute 
between coexisting gas and liquid, it is defined by 

kn = kBTn[/x = fcsTnfx exp(Az/s), (10) 

where kBTn[ is the solute partial pressure in gas and 
X = fiijn^^ is the solute molar fraction in liquid. In 
water in the ambient conditions, Ai/g is 3.44 for O 2 , 4.12 


I'g^n) = pg{n)/kBT. 
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FIG. 1: (Color online) Illustration of experimental setup with 
a pressure valve realizing a constant liquid pressure p (outside 
bubbles) larger than pcx, where the total solvent and solute 
numbers are Hxed. The cell contains bubbles in bulk (left) or 
those on a wall (right). The volume of the cell is V and that 
of the valve region is 14, where 14 is nearly equal to cpV for 
small bubble volume fraction tj) in the cell. 


for N 2 , and 0.18 for CO 2 , where CO 2 is highly soluble in 
liquid water. Thus, our theory is not applicable to CO 2 . 

However, there are a variety of solutes with stronger 
hydrophobicityjs^. For example, Ai/g = 10.2 for penta- 
cosane. In addition, from numerical simulations, a neu¬ 
tral hard-sphere particle deforms the surrounding hydro¬ 
gen bonding; as a result, Avg oc for a < 1 nm and 
Avg ^ /kfjT for a > 1 nm with varying the parti¬ 

cle radius a |2ll - l^ . This gives Ai^g ^ 180 for a ~ 1 nm. 
As hydrophobic assembly, such strongly hydrophobic so¬ 
lutes aggregate in liquid water. 


B. Chemical equilibrim and pressure balance 

We consider bubbles in bulk or on a wall at a small vol¬ 
ume fraction </> in the cell. For simplicity, we assume no 
bubble in the valve region in Fig.l. If the water density 
inside the bubbles n' is much smaller than h, the valve 
volume 14 is given by 


14 = V(j). (11) 

Since the total solvent and solute numbers are fixed, the 
densities in the liquid are given by 

h = n — (jin' = n, (12) 

ni = fli — (jin'i (13) 

Hereafter, we set n = n. We also have </> < ni/n{ ^ 1 
from fii > 0. Thus, the chemical equilibrium condition 
(8) and the conservation relation (13) give 

n'l = ni/[(j) + exp(-Azzs)], (14) 

h-i = hi/[l-I-(()exp(Ai/g)]. (15) 


The fraction of the solute in the bubbles is given by 

a = (jm'i/fii = (j)/[(j) + ex.p{-AKs)], (16) 

which tends to 1 for ^ ^ exp(—Ai^g). 

We write the value of the water chemical potential p 
in eq.(3) in the gas as p' and that in the liquid as fl, 
where p! = fivn equilibrium. Since p is close to Pcx, it is 
convenient to measure them from the chemical potential 
/icx on CX for pure water. Here, n' is small in the gas 
and use can be made of the Gibbs-Duhem relation in the 
liquid. Then, we obtain 

p' = Pc^ + kBT\\n{n'/nl^+n'^gs{n')], (17) 

p = 4 = /Xcx + (p-Pcx)/'«L> ( 18 ) 

where p remains equal to the initial value p. To linear 
order in the deviation n' — in eq.(17), the chemical 
equilibrium condition p' = p yields 

1 = {P - Pck)/ ksTni^ - n'^gsjn'), (19) 

In the right hand side of eq.(19), we may neglect the 
first term for p — pcx “C ksTn^^ and the second term for 
n'^vi <C 1 (see the sentence below eq.( 6 )). Then, we find 

n' = nl^ = Pcx/fesT. ( 20 ) 

For one-component fluids [s^, the pressure in a bubble 
has been set equal to Pcx from n^x/^^cx 1 from the 
critical point. In the present mixture case, the gas pres¬ 
sure is p' = ksTn'i +Pcx- With the aid of the Laplace law 
p' = p + 2(7/R, we obtain the pressure balance equation, 

kBTn'i=p-pcx + 2a/R. ( 21 ) 

Eliminating nj from eqs.(14) and (21), we may express 
the volume fraction (j) as 

^ = kBTnilip-p,^ + 2alR)-e-^''‘. ( 22 ) 

From eqs.(20) and (21), we find Uj » n' for p — pcx ^ 
Pcx or for R 2ajpcx, where the gas consists mostly of 
the solute. For water at T = 300 K, we have pcx = 3.6 
kPa and = 0.86 x 10^®/cm^, where n) 3> n' holds for 
p 0.0036 atm or for i? < 40 pvci. In addition, at p = 1 
atm, we have nj = 2a/kBTR for R ^ 1.4 pvci. 

In the limit of <()—>■ 0 and R — 00 , eq.(22) gives a 
threshold solute density for gas film formation, 

n'i = e-^''^{p-Pcx)/kBT, (23) 

which vanishes as p —> Pcx and is small for large Avg- 
Here, we introduce the following parameter, 

q = ni/ni-l. (24) 

A gas film can appear for 7 > 0, but bubbles with R~^ > 
0 can be stable for 7 > 7 tr with 7 tr being a positive 
threshold (see Fig.4). For O 2 in water at T = 300 K, 
we have n\ = 0.78 x 10^®(p —Pcx) cm“^ with pressures 
in atm. The corresponding oxygen mole fraction is 2.3 x 
10 “'’(p-pcx)- 
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C. Gas film at fixed pressure 

We consider a gas film on a hydrophobic wall, where 
there is no contact between the wall and the liquid phase 
Setting R~^ = 0 in eq.(22), we obtain </> for 7 > 0 as 

4> = (fii - nf)kBT/{p - pcx) 

= 7 exp(-A:/s). (25) 

In this case, we have hi = hi/(l + 7) from eq.(15). For 
7 < 0, water is only microscopically depleted at the wall, 
though the depletion layer itself can be influenced by the 
solute [13 ■ In the present isobaric case, </) increases and 
even approaches unity as p — Pcx ^ 0 , so we need to 
require p — Pcx > ksT^fii — nf). In contrast, at fixed 
cell volume, (j) remains small even for p — pcx A 0 (see 
Appendix A). 


D. Bubbles with a common radius at fixed pressure 


We suppose bubbles with a common curvature R~^ 
outside CX. Using the relation ksTfiiKp — pc^) = (7 + 
1 ) exp(—Ai/g), we rewrite eq.( 22 ) as 


kBT{ni — nf) R — Rc 

P-Pcyi R + iRc ’ 


(26) 


which tends to eq.(25) in the limit R 00 . Here, we 
introduce the critical radius Rc defined by 


i?c = 2 (t/7(p-Pcx)- (27) 

Here, Rc = l-d/y ^m for ambient water (300 K and 1 
atm). We need to require R > Rc outside CX since 
^ > 0. See Rc for O 2 in water in Fi g.5( a). For bubble 
nucleation in one-component fluids [35|-l37|, the critical ra¬ 
dius is given by Rc = 2(jj (pcx — P) with p < pcx- 

We assume N\y bubbles in the cell neglecting bubble 
coalescence. Then, we express 4> as 

(j) = 47ri?3G(6»)nb/3, (28) 

where nb = N\,/V is the bubble density. For bulk bub¬ 
bles we set G{9) = 1. For surface bubbles it is given by 
Turnbull’s formula [H, [s^l , 

G( 6 ») = (2-3 cos 6i-bcos3 6»)/4, (29) 

where 9 is the (gas-side) contact angle in the partial dry¬ 
ing condition determined by Young’s relation, 

cos 6 » = (crt - cr®)/cr. (30) 

where cr^ tr^ are the free energies per area between 
the wall and the liquid and gas phases, respectively, and 
we assume — cr® | < a. Here, 0 < d < 7 r /2 for a 
hydrophobic wall and -k< 9 < -k for a hydrophilic 
wall. As 0 —>■ 0, we have the complete drying condition 


cr^ — CTw = ^t T = Td on CX. As 0 —>■ TT, the bub¬ 
bles tend to be detached from the wall, resulting in bulk 
bubbles. Experimentally, 9 for surface bubbles has been 
observed in a range of 10-30° [3 • Note that eqs.(26) and 
(28) constitute a closed set of equations determining the 
equilibrium radius R for each given hi, 0, and rib. 


III. BUBBLE FREE ENERGY 

A. Derivation using grand potential density 

In the geometry in Fig.l with a pressure valve, we 
should derive the equilibrium conditions of bubbles from 
minimization of the Gibbs free energy written as 

G = F (tS + (V Vb}p 

= G + AG, (31) 

where F is the Helmholtz free energy (excluding the sur¬ 
face contribution here) and S is the total interface area. 
For a small volume change 14 —>■ 14 -I- dt4, the work 
exerted by the fluid to the valve is pdVh at fixed pres¬ 
sure, so we should consider G in eq.(31). The second 
line is the definition of the bubble free energy AG with 
G = Vf{n, hi) -I- Vp being the initial Gibbs free energy. 

In terms of the Helmholtz free energy densities f in 
the gas and / in the liquid, we have F = V {(j)f' + f) for 
the total system including the valve region. Here, it is 
convenient to introduce the grand potential density, 

uj(n,ni) = f - pn-pini+p, (32) 

where p and pi are the initial chemical potentials for 
water and solute, respectively. Using eqs.(12) and (13) 
we obtain 


F/V = + iP — (l)p + f, (33) 

where w' is the value of lo in the gas, P is that in the 
liquid, and / is the initial Helmholtz free energy density. 
Thus, eq.(31) gives 

AG = V[(j)uj' -I- A] + aS. (34) 

We note that w(n, ni) vanishes in the initial state and 
is second order with respect to the deviations n — n and 
ni—hi (see eqs.(12) and (13)). In the following we assume 
that the bubbles have a common curvature R~^, where 
in terms of G{9) in eq.(29) S is given by [ssL [^. 

S = A'KR^G{9)nbV. (35) 

We next calculate uj' and P for small p — Pcx assuming 
eqs.(14) and (20). In the gas, we use oj' = {fi' — p)n' + 
{pb'i — pi)n[ + p — p', where the first term in the right 
hand side is negligible from eq.(18). Further we set p' = 
Pcx + kBTn'j from eq.(20)and use eq.(3) to find 

ut' = kBTn[\b\{n[/ni) - 1 - Ai/^] + p - Pcx- (36) 




5 


oxygen = 0.149 \p.m '^] 



FIG. 2: (Color online) Normalized bubble free energy 
AG/N\^kBT vs R for O 2 in liquid water in the isobaric condi¬ 
tion for bulk bubbles in (a) and (b) and for surface bubbles in 
(c) and (d), where T = 300 K and rib = N\,/V — 0.149//rm®. 
The pressure difference Ap — p — pcx is 1 atm in (a), (c), 
and (d) and is 300 atm in (b). The contact angle 0 is 15° 
in (c) and is 6° in (d). Then, (A, 7m) = (3.86,7.74) in (a), 
(0.0129, 0.0680) in (b), (0.367, 0.945) in (c), and (0.109, 0.353) 
in (d), while 7tr — 1.147m for all the cases. In the equi¬ 
librium state at R = R 2 with AGmin < 0 on lowest curve 
in each panel, {■y,(fi,R 2 ) = (9.5,4.25 x 10“^, 0.41) in (a), 
(0.083,1.45 X 10"®, 0.13) in (b), (1.2,1.62 x 10■^3.1) in (c), 
and (0.44,7.00 x 10"®, 7.9) in (d) with R 2 in pm. 


In the liquid, hi is very small and we obtain 

a) = A:Br[hiln(hi/hi) -|- t/fnj], (37) 

where hi — hi — (j)n[ from eq.(13). Thus, if </) <C hijn^ <C 
1, the logarithm ln(fii/hi) can be expanded with respect 
to (/), leading to w oc However, we are also interested 
in the case (j) ~ hijn[. 

In equilibrium, AG in eq.(34) is minimized with re¬ 
spect to n[ and (j), where R and S are functions of (j) at 
fixed 0 and rib. Thus, let us change n[ and (p infinites¬ 
imally by Sn[ and 6(j), respectively. From eqs.(34)-(37), 
we calculate the incremental change of AG as 

5{AG) = VkBT[\ii{n'i/hi) — Ai/slin'iSp + (pSn'i) 

+V[p — Pcx + ‘2cr/R —kBTn{]6(j). (38) 


Therefore, the equilibrium conditions (14) and (21) follow 
from d{AG)/dn[ = d{AG)/d(p = 0. 

Furthermore, if we assume the pressure balance (21) 
(without assuming eq.(14)), AG becomes a function of R 
only under eqs.(28) and (35). Its derivative with respect 
to R is calculated as 


d{AG) 



- Av. 


4cr 


(39) 
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FIG. 3: (Golor online) Two radii Ri and R 2 vs 7/7m giving 
the local maximum and minimum of AG for O 2 in water at 
T = 300 K in the isobaric condition. They are written in 
dotted lines in the region 7m < 7 < 7tr. As in Fig. 2, displayed 
curves are for bulk bubbles with Ap = p — Pcx being 1 or 300 
atm (left) or for surface bubbles with 9 = 15° or 6° (right). 
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FIG. 4: (Color online) (a) 7 m(A) and ratio 7 tr(A)/ 7 m(A) vs 
A in the isobaric condition, where AG{R, 7, A) as a function 
of R exhibits local extrema for 7 > 7m and its local minimum 
becomes negative for 7 > 7tr. Here, 7tr/7m = 1.14 for any 
A. (b) Local maximum AGmax( 7 ,A) at R = Ri and local 
minimum AGmin (7, A) at R — R2 divided by AGm(7m, A) as 
functions of 7/7m for A = 1. Here, AGmin < 0 for 7 > 7tr, 
where R 2 is the equilibrium bubble radius. 


The extremum condition d{AG)/dR = 0 gives n[ = 
hie^''‘, leading to eqs.(14) and (15). 


B. Local maximum and minimum of bubble free 
energy at fixed pressure 

In Fig.2, we plot AG/N^kBT vs R for O 2 for bulk and 
surface bubbles in water under eq.(21), where T = 300 
K and p — Pcx = 1 or 300 atm. When we use O 2 (in 
Figs.2, 3, 6, and 7), we fix the bubble density at rib = 
Nh/V = 0.149//rm®. We recognize that AG assumes a 
local maximum at i? = i?i and a negative minimum at 
R = R2 for sufficiently large hi. Therefore, bubbles can 
appear in equilibrium at i? = i ?2 with increasing hi. For 
the case p — Pcx = 300 atm, the pressure in the bubble 
interior is also nearly equal to 300 atm for i? ;g> 5 nm 


dR 
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and the interior oxygen density is n[ = 7.2 /nm^ from 
the ideal-gas formula {p = nksT). Instead, if we use 
the van der Waals equation of state {p = nkBT/{l — 
n/3nc) — {9kBTc/8nc)n'^) at p = 300 atm and T = 300 
K, the density becomes 7.9 /nm^, where = 154.6 K 
and ric — 8.0 nm“^ for 02- Thus, the van der Waals 
interaction among O2 molecules is smaller than (per 
molecule) even at p — Pcx = 300 atm. 

To explain Fig.2, we treat AG as a function of R by in¬ 
creasing hi or 7 with the other parameters fixed. As will 
be shown in Appendix B, AG monotonically increases 
for 7 < 7m and exhibits a local maximum AGmax at 
R — Ri and a local minimum AGmin at R = R 2 , where 
AGmin > 0 for 7in < 7 < 7tr and AGmin < 0 for 7 > 7tr. 
Here, R 2 /R 1 increases from 1 with increasing 7 above 
7m. In each panel in Fig.2, we set 7m < 7 < 7tr for the 
upper curve, 7 = 7tr for the middle curve, and 7 > 7tr 
for the lower curve. Therefore, the two-phase states at 
R = R 2 are metastable for 7m < 7 < 7tr and stable for 
7 > 7tr- In Fig.3, we plot Ri and R 2 for O2 in water for 
bulk and surface bubbles. 

In Appendix B, we shall see that 7m and 7tr depend 
only on the following dimensionless parameter, 

A = [2cr/(p-pcx)][47rG(6>)nbe^''V3]^^^, (40) 

which diverges as p ^ Pcx and becomes small with in¬ 
creasing p — Pcx and/or decreasing G{9)n\^. Using A, we 
may rewrite eq.(26) in terms of u = R/Rc as 

= ^^{u-l)/{u + j), (41) 

which holds for u = Ri/Rc and i?2/7?c- In Fig.4, we plot 
7m and 7tr vs A in (a) and display AGmax and AGmin as 
functions of 7/7m at A = 1 in (b). Here, 7m ~ A^/^ for 
A <C I and 7m ~ A for A I. For any A, we find 

7tr/7m = 1A4. (42) 

For example, 7m = 0.330,2.23, and 19.4 for A = 0.1, 1, 
and 10, respectively. The threshold of bubble formation 
(7 > 7tr) is thus approximately given by hi > n\ for 
A < 1 and by hi > Arij for A > 1. 

In particular, with increasing the solute density, we 
examine the case 7 ^ 7m ~ 7tr using eq.(41), where R 2 
is the equilibrium radius. Then, for any A, we find 

7^1 + A3(1+7)774+ ...^ 1. (43) 

The ratio Ri/Rc rapidly approaches 1 with increasing 
7. In fact, even at 7 = 7tr, we have Ri/Rc = 1.091, 
1.122, and 1.150 for A = 0.1, 1, and 10, respectively. On 
the other hand, supposing u = R 2 /RC ^ 1, we obtain 
u^{u -I- 7) = 74/A^ from eq.(41). For A ^ 1, we have 

R 2 /RC ^ > 1. (44) 

For A ^ 1, there are two limiting cases: 

R 2 /RC ^ (7/A)3/2>1 (A<7<A3) (45) 

^ 74/VA>7 (7>A3). (46) 


In these limiting cases, we surely obtain R 2 ^ Rc- See 
Fig. 10 in Appendix B for the behaviors of i?i and i?2 vs 
7. On the other hand, for 7 ^ 7m) the solute fraction a 
in bubbles in eq.(16) is much smaller than 1 at i? = i?i 
and approaches 1 at i? = i?2 (see Fig. 10(b)). 

C. Bubble free energy at fixed chemical potentials 

So far we have fixed the total particle numbers N = 
V{h + (j)n') and Ni = V{hi + (j)n'i) as well as the liquid 
pressure. In this case, the water chemical potential jx is 
nearly fixed at the initial value fl from eq.(18). As an¬ 
other boundary condition, we may attach a solute reser¬ 
voir to the cell to fix the solute chemical potential at the 
initial value pi, where we still attach a pressure valve. In 
this grand canonical case, we have nj = hi so that A = 0 
from eq.(32). We should minimize the grand potential, 

H = G — fiN — p-iNi 

= V(I)uj' + <tS, (47) 

where G is defined in eq.(31), p, and pi are the initial 
chemical potentials, and uj' is given by eq.(36). To de¬ 
rive the second line of eq.(47), we have used the relation 
F/V = (/w' — (1 -b (l))p + p{n + (j)n') + pi{ni + (pn[), where 
N and Ni are not fixed. 

With respect to small changes —>• -b Sn'i and (j) 

4> -b <50, the incremental change of H is given by the right 
hand side of eq.(38) if hi is replaced by hi. Therefore, 
the extremum conditions dVl/dn'i = dfl/dcj) = 0 yield 
the pressure balance (21) and the chemical equilibrium 
condition n[ = hie^"^®. If these extremum conditions are 
assumed, we obtain 

w'=-(p-Pcx)7) (48) 

which is negative for 7 > 0 outside CX. Here, in the 
second line of eq.(47), the first term is proportional to 
R^ and the second term to R^, so the minimum of H 
decreases monotonically with increasing R (for R> Rc), 
indicating appearance of macroscopic bubbles. 

IV. SOLUTE-INDUCED NUCLEATION 
A. Experimental situations 

We have shown that the bubble free energy AG{R) has 
a local maximum at R = Ri and a minimum at R = R 2 
for 7 > 7tr (except for gas films). In such situations, the 
initial homogeneous state is metastable and there can 
be solute-induced bubble nucleation outside CX (s^]. In 
contrast, in one-component fluids, bubble nucleation oc¬ 
curs only inside CX {p < pcx) |35l439l| . In nucleation, cru¬ 
cial is the free energy Fc needed to create a critical bub¬ 
ble with R = Rc- We call it the nucleation barrier, since 
the nucleation rate / of bubble formation is proportional 
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FIG. 5: (Color online) (a) Critical radius Rc in eq.(27) vs 
oxygen density ni for bulk bubbles in water, where Ap = p — 
Pcx = 0 , 10 , and 300 atom and T = 300 K. (b) AGone{R)/Fc 
in eq.(51) vs R/Rc for 7 = 0.5,4, and 49. 


oxygen rib = 0.149 [^m 



FIG. 6 : (Color online) Normalized nucleation barrier Fc/ksT 
vs 7 for O 2 in water, for (a) bulk and (b) surface bubbles. 


to the Boltzmann factor exp(—Fc/fc^T). Therefore, if 
Fc/ksT is too large (say, 80), I becomes too small for 
experiments on realistic timescales. In our case, is 
reduced with increasing p — pcx and/or 7 = hi/nf — 1 . 
For surface bubbles, it is also reduced with decreasing 
the contact angle 9. 

We make some comments on experimental situations. 
First, in the previous observations (isl - ll^ . bulk nanobub¬ 
bles have been produced by breakup of large bubbles 
composed of a gas such as O 2 , CH 4 , or Ar, where the 
typical flow-induced bubble size is of great interest [4l| . 
Second, a small amount of surfactants and/or ions are 
usually present in water, which increase the bubble sta¬ 
bility 0. Indeed, surfactant molecules at the gas-liquid 
interface reduce the surface tension, while electric charges 
or electric double layers at the interface prevent bubble 
coalescence [dol - li^ . For example, with addition of O 2 
and a salt in water [l^. the bubble-size distribution on 
long timescales was found to have a peak at i? ^ 100 
nm. Third, on a non-smooth hydrophobic wall, there can 
be preexisting trapped bubbles or strongly hydrophobic 
spots. In such cases, there should be no significant nu¬ 
cleation barrier for the formation of surface bubbles with 
small contact angles 9. 


B. Critical radius and nucleation barrier 


We consider a single bubble with curvature R~^ in 
bulk or on a hydrophobic wall. In the early stage with 
small 4>, we may neglect w cx in eq.(34) to obtain the 
single-bubble free energy in the standard form [ssl - fs^ . 


AGo„e(i?) = G{9) 


—ttR^u}' + AiraR^ 
3 


(49) 


is maximized at the critical radius Rc in eq.(27). See 
Fig.5(a) for Rc vs hi for O 2 . Since w' = —^ajRc at 
R — Rc from eq.(36), the nucleation barrier (=the max¬ 
imum of AGone at i? = Rc) is written as 

Fc = 47rG(6>)cri?//3 = 167rG(6»)cr^/[(p - Pcx)7]^- (50) 

For surface bubbles with small d, we have G(9) = 30^/16 
and Gc oc For bubble nucleation in one-component 
fluids, Fc is given by the above form with 7 = 1 and 
P < Pcx- In terms of u = R/Rc, we may also express 
AFone{R) simply as 

AGone/Ac = 2u^{uIj + 1 ) In (51) 

The right hand side may be approximated by —2u^ + 
‘iv? for 7 <C 1 and by —2u^lnu -I- for 7 ^ 1. In 
Fig.5(b), we plot the above scaling function. In addition, 
the nucleation rate / is of the form [s^ 

I = Tcni,eM-Fc/kBT), (52) 

where Fc is the growth rate of a critical bubble (see 

eq.(60) in the next subsection). 

For water at T = 300 K, we have 

Fc/kBT^73G{9)Rl (53) 

with R in nm. In Fig. 6 , we plot Fc/kBT vs 7 for O 2 in 
water for bulk and surface bubbles. . In homogeneous 
bubble nucleation of pure water at T ~ 300 K , 

bubbles with R > Rc are detectable for Fc/kBT < 70 or 
for Rc 7 nm in experimental times and Rc can be of 
order 1 nm only for negative p of order —1000 atom. For 
O 2 in our case, Rc is decreased down to 1 nm, depending 
on hi, p — pcx, and 9 in Fig. 6 . 


where w' is given by eq.(36) and nj is related to R by 
the pressure balance (21). Note that w' is usually a 
negative constant in nucleation in metastable systems. 
Here, d{AGc,ne)/dR cx ln(n{/hi) — Az/g, which follows 
from eq.(39) if hi is replaced by hi. Then, AGone(7?) 


C. Dynamics of bulk nucleation 

We next examine nucleation dynamics of bulk bub¬ 
bles for 7 > 7 tr. To describe attainment of the equi¬ 
librium radius i ?2 in the simplest manner, we assume a 











Bubble growth in bulk, oxygen, nb = 0.149 [/um 




FIG. 7: (Color online) Growth of radius R{t) for bulk bubbles 
in the isobaric condition for O 2 in water at T = 300 K, which 
are obtained from eq.(59) for R{0) — Ri + 0.3 A with D = 
2 X 10“® cm^/s. In (a) and (a’), 7 = 10 and 15 with Ap = 
p—Pcx =1 atom. In (b) and (b’), 7 = 2 and 7.5 with Ap =300 
atm. Curves are written on linear scales in (a) and (b) and 
on semi-logarithmic scales in (a’) and (b’). Here, 7 tr = 8.80 
and 0.0774 for Ap = 1 and 300 atm, respectively. 


common radius R{t) for all the bubbles with a constant 
bubble density rib. We also assume a time-dependent 
background solute density in the liquid defined by 

hi(t) =ni- (54) 

where (j){t) is determined by eq.(28) with G{9) = 1. Ex¬ 
pressing (j)(t) and nj(t) in terms of R{t), we may describe 
saturation of (/)(t) up to the equilibrium volume fraction. 
After this stage, however, the bubble number decreases 
in time in the presence of bubble coalescence (which can 
be suppressed with addition of salt [dol. l4ll| ). 

For simplicity, we further assume that the solute dif¬ 
fusion constant D is much smaller than the thermal dif¬ 
fusion constant Dp in the liquid. Then, we can neglect 
temperature inhomogeneity around bubbles, which much 
simplifies the calculation. In fact, for liquid water at 300 
K and 1 atom, we have Dp 1.4 x 10“^ cm^/s and 
D - 2.0 X lO-'^ cmVs (< Dp) for O 2 . 

We focus our attention to a single bubble neglecting its 
Brownian motion, where nj (r, t) slowly changes in time t 
tending to ni{t) in eq.(54) far from it. We write the dis¬ 
tance from the droplet center as r. In the bubble exterior 
r > R, the solute obeys the diffusion equation. 
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We assume the continuity of the solute chemical potential 
fii aX r = R + 0 and r = R — Q across the interface. From 
eq.(7), the solute density nf = ni{R + 0,t) immediately 
outside the bubble is related to the interior density n[ by 

(56) 

Therefore, in the quasi-static approximation ( 3 ^ . ni(r,t) 
slightly outside the interface is written as 

ni(r, t) = fii + (np — ni)R/r. (57) 

The flux to the bubble is given by D(ni — np)/R, so the 
conservation of the solute yields 

(n[ -np)^ = ^(ni- np) (58) 

Here, ni — nf = fii— n((r/>+e~^''^) from eqs.(13) and (55, 
so the right hand side of eq.(58) vanishes at R = R 2 from 
eq.(13). In accord with the equilibrium relation (27), 
division of eq.(58) by n[D/R gives the desired equation, 




R dR 

’15~di 


(ni - ni)(l - Rc/R) 
(p -Pcx + 2alR)lkBT 


</). (59) 


For 7 > 7 tr, the right hand side of eq.(59) vanishes for 
R = Ri and R 2 , where Ri = Re- Here, bubbles with 
R > Ri grow up to i? 2 , while those with R < Ri shrink. 
If the deviation 6R = R — Ri is small, it obeys the linear 
equation d(6R)/dt = T^SR, where T^ is the growth rate 
of a critical bubble of the form, 

r, = Z4i?-V[(l + 7)(e^"^-l)]- (60) 


In terms of and u = R/Rc, we may rewrite eq.(59) as 


— = Tc- - - 

dt u \_u + j 7 ^ 


(61) 


which is consistent with eq.(41). 

In Fig.7, we display the growth of R(t) by setting D = 
2 X 10“® cm^/s for O 2 in water at T = 300 K, where 
P — Pcx is 1 atm in (a) and (a’) and 300 atm in (b) and 
(b’). As the initial radius, we set i?( 0 ) = Ri + 0.3 A, 
which yields u(0) — Ri/Rc ^ 10“^ in eq.(61). The right 
panels indicate the exponential growth, 

i?(t) = i?( 0 )-b(i?( 0 )-i?i)e^'=‘ (62) 


in the early stage. Numerically, Tc is 2.89 and 6.70 for 
7 = 10 and 15, respectively, in (a) and (a’), while it is 
7.63 and 142 for 7 = 2 and 7.5, respectively, in (b) and 
(b’). These values agree with eq.(59). In this calculation, 
we assume the pre-existence of bubbles with radii slightly 
exceeding Ri. However, if we start with the homogeneous 
initial state, the birth of such large bubbles in the cell 
occurs as rare thermal activations on a timescale of order. 


(55) 1/VI - exp(F,/fcBT)/F<,r,. 


dni 


DV^m. 


(63) 
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V. SUMMARY 

We have investigated bubble formation in bulk and on 
hydrophobic walls induced by accumulation of a small 
amount of a neutral solute in liquid water outside the 
solvent CX. We have used the fact that a gas such as O 2 
or N 2 remains in gaseous states within phase-separated 
domains in ambient liquid water, because it is mildly hy¬ 
drophobic with a critical temperature much below 300 
K. With this input, we have constructed a simple ther¬ 
modynamic theory for dilute binary mixtures including 
a considerably large solvation chemical potential differ¬ 
ence A^s = ksTAvs- We have assumed fixed particle 
numbers and a fixed liquid pressure (N-Ni-p) in the text 
and in Appendix B, but we have also treated bubble for¬ 
mation in the fi-pi-p ensemble in Sec.IIIC and in the 
N-Ni-V ensemble in Appendices A and B, 

In particular, in Sec.II, we have found a threshold so¬ 
lute density rij in eq.(23) for film formation on a com¬ 
pletely dried wall at fixed pressure, which is very small 
for large At'g. The threshold density is increased to 
(7m + l)?^i for metastable bubbles and (ytr + l)nf for 
stable bubbles due to the surface tension, where ytr — 
1.147 i„. Here, 7 i„ and 7 tr are displayed in Fig.4(a) as 
functions of a parameter A in eq.(40). In Sec.Ill, we 
have also presented a bubble free energy AG for a small 
gas fraction (f) in eqs.(34)-(37) for the isobaric case, whose 
minimization yields the equilibrium conditions (14) and 
(21). In Sec,IV, we have calculated the critical radius 
Rc and the barrier free energy for solute-induced nu- 
cleation. The Fc/ksT is very high for homogeneous nu- 
cleation except for high liquid pressures, but it can be 
decreased for heterogeneous nucleation with a small con¬ 
tact angle 6. 

We make some critical remarks, (i) First, we have as¬ 
sumed gaseous domains. However, with increasing p—pcx 
and/or Aug, liquid or solid precipitates should be formed 
in bulk and on walls, sensitively depending on their mu¬ 
tual attractive interaction. Note that a large attractive 
interaction arises even among large hard-sphere particles 
in water due to deformations of the hydrogen bonding 
[2ll - [^ . (ii) Second, we should include the effects of sur¬ 
factants and ions in the discussion of the bubble size dis¬ 
tribution In this paper, we have assumed a 

constant bubble density rib in eqs.(28), (35), and (40). 
This assumption can be justified only when bubble co¬ 
alescence is suppressed by the electrostatic interaction 
near the gas-liquid interfaces, (hi) Third, dynamics of 
bubble formation and dissolution!^ should be studied 
in future, which can be induced by a change in pressure, 
temperature, or solute density. 
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Appendix A: Bubbles at fixed cell volume 

Here, we consider two-phase coexistence outside CX, 
fixing the particle numbers and the cell volume V without 
a pressure valve. Some discussions were already made 
on attainment of two-phase equilibrium in finite systems 
inside CX (without impurities) [13, ■ 

In this case, while eqs.(14) and (15) are unchanged, 
the liquid volume is decreased by 4>V for n' <C h with 
appearance of a gas region. As a result, the liquid density 
fi is increased as 


n = (1 -I- 4>)n. 


(Al) 


In terms of the isothermal compressibility Kt of liquid 
water, the pressure increase is given by so the 

pressure balance relation ( 21 ) is changed as 


ksTn'i = p - pcx +(I>/Kt+ ‘2(t/R. (A2) 

Here, Kt = 0.45 x 10“^/MPa in ambient water near CX, 
where even a very small (j) gives rise to a large pressure. 

The bubble free energy AF is defined as the increase 
in the Helmholtz free energy as F = F -|- AF due to 
appearance of bubbles. Some calculations give 

AF = V[(l)uj' + {1- ct>)uj] + aS. (A3) 

Here, we may replace (1 — (f>)6j by uj for small (j). Then 
AF assumes the form of eq.(34), but we need to change 
w in eq.(37) as 

a) = T[ni ln(ni/fii) -|- (j)n[] + /2Kt, (A4) 


where the last term is due to the compression in the liq¬ 
uid. From eq.(l), it is equal to /"(n)(n — n)^/2, where 
/" = = \Iv?Kt- The counterpart of eq.(38) 

for the increment d{AF) is obtained iip — pcx is replaced 
by p — Pcx + Minimization of AF with respect to 

(j) and n[ thus yields eqs.(14) and (A2). Furthermore, the 
derivative d{AF) /dR is obtained if p —pcx is replaced by 
P — Pcx + 2(1)/Kt in the right hand side of eq.(39). 

The equilibrium equation for i? or is given by eq.(22) 
if p — Pcx is replaced by p — pcx + 4>/Kt- In particular, 
for a gas film (i?“^ = 0 ), (j) is explicitly calculated as 


= e 


— Aus 


\/(l J- /i)2/4 -I- 7/1 


1 + h 
2 


(A5) 


where we define 


h = Kt{p- Pcxje^''^- (A 6 ) 

For O 2 in ambient water, we have h = 1.4 x 10“^(p —Pcx) 
with pressures in atm, so <C 1 for p ^ 10 ^ atm. 

Here, we assume 1 J- /i > 0. A film appears for 7 > 0 as 
in the isobaric case. In particular, for 7 |/i| <C (1 J- ^)^/4, 
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we find the linear behavior oc 7 as 


= 76 


— Aus 


Kx 
1 + h 


h/{l + h) 

ksTnie^''^ -p + Pcx 


(A7) 


From the first line, this formula tends to eq.(25) only for 
h ^ 1. The second line can be used even for p < pcx, 
where 4> increases with increasing pcx ~ P and/or hi. See 
Appendix B for more analysis for the case R~^ > 0. 

Appendix B: Scaling of bubble free energy 


Here, we examine the bubble free energy AF in eq.(34) 
by scaling it in a dimensionless form, assuming a common 
curvature R~^ for all the bubbles. 

1. Fixed pressure 

At fixed pressure in Fig.l, we assume the pressure bal¬ 
ance ( 21 ) and introduce scaling variables s and v by 

s = (j)e^''‘ = (47rG(6»)nb/3)e^‘^=i^^ (Bl) 

V = [2a/(p-pc^)R]= As~^/^. (B2) 

where A is the parameter in eq.(40). From eq.(16) the 
solute fraction in bubbles is a = s/(l + s). As a scaled 
bubble free energy, we define F as 

F=e^''‘AG/[V{p-p,^)]. (B3) 

From eqs.(34)-(37), we express F in terms of s and v as 


= [7 -b 1 — s(l -b u)] In 
1 -b u 


1 — s 


1 -b u 


1 + 7. 


-bs(l -b v) In 


.1 + 7. 


-b + 2 ®'^ + 


(B4) 


where 7 is given by eq.(24). With fixed 7 and A, F is 
a function of s only. From eq.(B2) its derivative with 
respect to s is calculated as 


+.F=(l + |„|l„ 


1 -b u 


.7 + 1 - s(l -b v) 
The extremum condition dF/ds = 0 yields 
7 = s -b (1 -b s)v = s -b A (1 -b 


(B5) 


(B 6 ) 


which is equivalent to eqs.(21), (26), and (41). If eq.(B 6 ) 
is assumed, we have F = (3u/2 -b l)s — (1 -b 7 ) ln(l -b s) 
as extremum values depending only on A. In Fig. 8 , we 
display F{j, s, A) in the 7 -s plane at A = 1. 

For each A, the right hand side of eq.(B 6 ) is minimized 
at s = Sm(A) as a function of s, where Sm satisfies 

A = 34/V(1 - 2sm). (B7) 



FIG. 8 : (Color online) F in the isobaric condition in the 7 -s 
plane at A = 1 in the neighborhood of 7 = 7m and s = Sm. 
Curve of {dF/ds)-y = 0 or eq.(B 6 ) is written on the surface, 
on which s = si (7, A) for the local maximum and s = S2 ( 7 , A) 
for the local minimum for each 7. 



FIG. 9: (Color online) Results in the isobaric condition, (a) 
7 / 7 m vs s/sm from the extremum condition (B 6 ) for A = 0.1, 
1, and 10 and (b) Sm and Str vs A. Here, 7 m, Sm, and Str are 
determined by eqs.(B7)-(B9). 


As can be seen in Fig.9(a), if 7 > 7 „i, eq.(B 6 ) has 
two solutions 31 ( 7 , A) and 32 ( 7 , A) with 3 i < 32 , where 
.F(s , 7 , A) exhibits a local maximum Jqnax(7) A) at 3 = 3i 
and a local minimum Jqnin(7)A) at 3 = 32. Further in¬ 
creasing 7 above 7 ni, the local minimum Jqnin decreases 
and becomes negative for 7 > 7 ti.(A), where 7 tr and the 
corresponding 3 , written as 3tr(A), are calculated from 

A/(A + sli^) = 2(1 + 3 ,-/) ln(l + 3tr) - 2, (B9) 

7tr = Str + A(1 -b 3tr)/St/^- (BIO) 

See Fig.4(a) for 7 ^ and 7tr/7ni(— 1.14) vs A. 

From eqs.(B7)-(B10), we seek the asymptotic behav¬ 
iors for small and large A. For A <C 1 we find 

Srrr = {A/2,f/\ 3tr = A^/^ 

7m = 4(A/3)3/4, 7t, ^ 2A3/4. (Bll) 


The minimum of eq.(B6) at 3 = 3m is written as 


On the other hand, for A 1, we have 


(B 8 ) 3m ^ 1/2, 3tr = +401 


7m — Sm(4 -b 3m)/(1 23m). 
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7/7m 7/7m 



sure balance (A2) and use w in eq.(A4). Introducing the 
scaled bubble free energy T as in eq.(B3) (with replace¬ 
ment AG —> AA), we express it in terms of s in eq.(Bl), 
V in eq.(B2), and 


V = {2a /R + (l)fKT)l{p - Pcx) =v + s/h, (B14) 
where h is defined in eq.(A6). Replacing u by D in the 
first two terms in eq.(B4), we obtain 


J- = ['y + 1 — s(l -I- u)] In 
1 -I- fi 


1 — s 


1 


- V 


-|-s(l -I- v) In 


.1 + 7. 


1 -f 7_ 
,2 


3 s* 


(B15) 


where the last term arises from the compression term in 
eq.(A4). As in eq.(B5), the derivative of F with respect 
to s is calculated as 


ds 


F= 1 


2v 2s 

7i" 


In 


.7+1 


The extremum condition dFjds = 


1 + v 

- s(l -I- v) 
0 yields 


(B16) 


FIG. 10: (Color online) Results in the isobaric condition, (a) 
Ri/Rm = (si/sm)+® and i? 2 /Am = (S 2 /Sm)+® VS "f /and 
(b) oi = si/(l-|-si) and 0:2 = S 2 /(l-|-S 2 ) vs 7/7m, where A = 
0.1,1, and 10. (c) Ri/Rc — 1 = sJ+y/A — 1 vs y/ym, which 
is small even for vs y = ytr (see eq.(43)). (d) i? 2 /[Rcy"^+/A] 
for A = 0.1, 1, and 10 (bold lines) and i? 2 /[Rc(y/A)®+] for 
A = 10 (dotted line) as functions of y/ym (see eqs.(44)-(46)). 

y^ ^ (3/22/3)A, ytr = 2.145A, (B12) 

where Str and ytr are calculated numerically. If y 3> ym, 
we obtain eqs.(43)-(46). 

For y > yni, we consider the radii Rm, 4?i, and R 2 
corresponding to Sin,si, and S 2 - From eq.(B2) we have 


y = s -b (1 -b s)v 

= s -b A(1 -b s)s“+3 _|_ (g _|_ (B17) 

As in the fixed pressure case, F{s,'y,A,h) exhibits a 
local maximum at s = si(y,A, h) and a local minimum 
at s = S 2 (y,A, h) for y > j^{A,h) and its local mini¬ 
mum becomes negative for y > ytr (A, h). Note that the 
right hand side of eq.(B17) is minimized at s = Sm(A, h), 
where Sm is determined by 

A = 3s+3[I + (1 + 2s^)/h]/{l - 2s^). (B18) 

The corresponding minimum of eq.(B17) is written as 
ym = [4 -b Sm -b 4(1 -b Sni)^/h]Sni/(1 — 2Srn) (B19) 


Am _ 1/3 Ai _ 1/3 R 2 

Ito rCo rto 


1/3 


(B13) 


where Rq = 2cr/[(p —pcx)A] = Rc'^/A. In Fig.10, we plot 
Al/Am = (si/Sm)^/^ and R2/Rm = (S2/Sm)^/^ VS y/ym 
in (a) and ai = si/(l-bsi) and a 2 = 52/(1 + 52 ) vs y/ym 
in (b), where the latter are the solute fractions in bubbles 
dX R = Ri and ^ 2 - Furthermore, in (c), Ri/Rc — 1 is 
shown to be small for y slightly larger ym in accord with 
eq.(42). In (d), we divide R 2 IRC by its asymptotic forms 
for y 3> ym to confirm eqs.(44)-(46). 

2. Fixed volume 

We next scale the bubble free energy AF in the fixed- 
volume condition in Appendix A. We assume the pres- 


Here, we assunme Ah <C 1, where 5m (Ah)^^^ <C 1 
from eq.(B18). In fact, for O 2 in ambient water, we ob¬ 
tain Ah = O.OlO[nbG(0)]3/3 with Ub in units of 
which is independent of p — Pcx- Then, as in eq.(Bll), 
we find the asymptotic behaviors, 

5m = (Ah/3)3/4, Str = {Ahf^'^ 

ym = 4(A/3)3/Vh+4, ytr ^ 2A^l^lh^l\ (B20) 

where ytr/ym — 1+4. As in eqs.(43) and (44), the radii 
Ai and R 2 behave for y ym as 

Ai ^ Am A 2 = lA > i?m (B21) 
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